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ABSTRACT 

We study the intermittency and field-line structure of the MHD turbulence in plasmas with very large magnetic 
Prandtl numbers. In this regime, which is realized in the interstellar medium, some accretion disks, protogalaxies, 
galaxy-cluster gas, early Universe, etc., magnetic fluctuations can be excited at scales below the viscous cutoff. 
The salient feature of the resulting small-scale magnetic turbulence is the folded structure of the fields. It is 
characterized by very rapid transverse spatial oscillation of the field direction, while the field lines remain largely 
unbent up to the scale of the flow. Quantitatively, the fluctuation level and the field-line geometry can be studied in 
terms of the statistics of the field strength and of the field-line curvature. In the kinematic limit, the distribution of 
the field strength is an expanding lognormal, while that of the field-line curvature K is stationary and has a power 
tail ~ K'^l 1 . The field strength and curvature are anticorrelated, i.e. the growing fields are mostly flat, while the 
sharply curved fields remain relatively weak. The field, therefore, settles into a reduced-tension state. Numerical 
simulations demonstrate three essential features of the nonlinear regime. First, the total magnetic energy is equal 
to the total kinetic energy. Second, the intermittency is partially suppressed compared to the kinematic case, as the 
fields become more volume-filling and their distribution develops an exponential tail. Third, the folding structure 
of the field is unchanged from the kinematic case: the anticorrelation between the field strength and the curvature 
persists and the distribution of the latter retains the same power tail. We propose a model of back reaction based 
on the folding picture that reproduces all of the above numerical results. 

Subject headings: galaxies: magnetic fields — ISM: magnetic fields — magnetic fields — MHD — turbulence 



1. INTRODUCTION 

In turbulent MHD systems where the ratio of fluid viscosity 
and magnetic diffusivity (the magnetic Prandtl number, Pr m = 
v/rj) is very large, there exists a broad range of sub viscous 
scales available to magnetic fluctuations, but not to hydrody- 
namic ones. This MHD regime is encountered, for example, in 
such astrophysical environments as the interstellar medium and 
protogalactic plasmas, where Pr m can be as large as 10 14 to 10 22 
(Kulsrud 1999). Since the ratio of the resistive and viscous cut- 
off wave numbers is k v /k v ~ PrJ( 2 , this gives rise to subviscous 
scale ranges 7 to 1 1 decades wide. 

Since the fluid is highly conducting, the magnetic-field lines 
are (nearly) perfectly frozen into the fluid flow. The fluid mo- 
tions, even though restricted to the scales above the viscous cut- 
off, can excite magnetic fluctuations at much smaller scales via 
stretching and folding of the field lines. This possibility was 
first indicated by Batchelor (1950). The weak-field (kinematic) 
limit has been an attractive object of analytical study since the 
seminal work of Kazantsev (1967). The spectral theory of the 
kinematic dynamo driven by a random velocity field predicts 
exponential growth of the magnetic energy and its accumula- 
tion at the resistive scales (see Kazantsev 1967; Kulsrud & An- 
derson 1992; Gruzinov, Cowley, & Sudan 1996; Schekochihin, 
Boldyrev, & Kulsrud 2002a, and references therein). More re- 
cently, it was realized that small-scale magnetic fields generated 
by this "stretch-and-fold" dynamo possess a distinctive spatial 
folding structure (Fig. 1): the smallness of the field scale is 
due to rapid transverse spatial oscillation of the field direction, 
while the field lines remain largely unbent up to the scale of the 



flow (Ott 1998; Kinney et al. 2000; Schekochihin et al. 2002c, 
— this last paper is henceforth referred to as SCMM02). 

With the dramatic increase in the reach of the numerical ex- 
periment, the nonlinear regime became increasingly amenable 
to detailed study. The pioneering work of Meneguzzi, Frisch, 
& Pouquet (1981) was in recent years followed by a number 
of numerical investigations (Cattaneo, Hughes, & Kim 1996; 
Brandenburg et al. 1996; Zienicke, Politano, & Pouquet 1998; 
Kinney et al. 2000; Cho & Vishniac 2000; Brandenburg 2001; 
Chou 2001; Brummell, Cattaneo, & Tobias 2001; Maron, Cow- 
ley, & McWilliams 2002, the latter paper henceforth abbrevi- 
ated MCM02). However, the physical difference between the 
Pr m = 1 and Pr m ^> 1 regimes is not always realized. 

Another defining physical feature of the MHD regime we are 
considering is the absence of an externally imposed uniform 
magnetic field. The difference is essential. First, a fixed uni- 
form field implies a nonzero net flux through the system. Sec- 
ond, in the presence of a strong such field, magnetic-field lines 
cannot be bent, so the physics of the subviscous-scale magnetic 
fluctuations is more akin to that of the scalar turbulence (cf. 
Cho, Lazarian, & Vishniac 2002). 

In the astrophysical context, the main question has been of 
the impact small-scale magnetic fluctuations have on the fea- 
sibility of generating the large-scale galactic magnetic field by 
means of turbulent dynamo. In particular, one wonders how the 
accumulated small-scale magnetic energy affects the applica- 
bility of the mean-field ail-dynamo theory, which, in one form 
or another, has been at the center of all attempts to build a the- 
oretical understanding of the galactic magnetic fields (Beck et 
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al. 1996). This issue remains unresolved in both the theoreti- 
cal and numerical senses of the word. We note that, because of 
the large scale separations that have to be captured, the brute- 
force numerical solution of this problem remains beyond the 
capacity of currently available computational resources. On the 
other hand, the physics of the small-scale magnetic turbulence 
can be effectively studied as a separate problem and accessed 
(if only just) by numerical experiment. We hope that by de- 
veloping a thorough physical understanding of the small-scale 
magnetic fluctuations, we can approach the problem of their in- 
terplay with the large-scale fields and motions. 

In this work, we concentrate on the structural properties of 
the small-scale magnetic fields: namely, the intermittency of 
their spatial distribution and the geometry of the field lines. 
Quantitatively, these are studied in terms of the one-point statis- 
tics of the field strength and of the field-line curvature. 

We consider the nonlinear turbulent dynamo to be described 
by the equations of incompressible MHD: 
d 

— u = ^Au-V/? + B-VB + f, (1) 

at 

d 

— B = B Vu + r/AB, (2) 
dt 

where d/dt = d, +u • V is the convective derivative, u(f ,x) is the 
velocity field, B(f ,x) is the magnetic field, and f(f,x) is a ran- 
dom (white in time) large-scale forcing. The density p of the 
plasma is taken to be constant. The incompressibility condition 
V • u = is, therefore, added to the equations above and serves 
to determine (or, indeed, to define) pressure. For the sake of 
convenience, the pressure p and the magnetic field B have been 
normalized to p and (4-Trp) 1 / 2 , respectively. 

The plan of further proceedings is as follows. In Sec. 2, we 
review the necessary facts about the kinematic regime of the 
dynamo. These are important, for they form the basis of our 
understanding of the dynamo and remain surprisingly relevant 
in the nonlinear regime. The latter constitutes our main object 
of study. In Sec. 3, we describe the results of our numerical ex- 
periments and propose a heuristic physical model that explains 
these results. Conclusions are drawn in Sec. 4. 




FIG. 1. — Folding-structure formation via stretching of the field lines. Bold 
arrows indicate directions in which volumes are stretched by random shear. 

2. THE KINEMATIC REGIME 

In the weak-field (kinematic) limit, MHD turbulence reduces 
to the problem of passive advection of a vector field by a tur- 
bulent velocity field. The magnetic energy grows exponentially 



and the small-scale folding structure is formed at the time scale 
associated with the eddies that turn over the fastest (i.e., in Kol- 
mogorov turbulence, the viscous-scale eddies). Physically, this 
follows from the fact that the turbulent eddies act on the small- 
scale fields as a sequence of random linear-shear transforma- 
tions (see Fig. 1 and §2.3). An expanding lognormal distri- 
bution of the field strength emerges, which is qualitatively ex- 
plained in terms of the Central Limit Theorem. 

2.1. The Kazantsev-Kraichnan Model 

All of the above results can be derived analytically in the 
framework of the Kraichnan (1968) model of passive advection, 
which replaces the turbulent velocity with a Gaussian random 
field (^-correlated in time: 

( M ''(f,xV(f',x')) = <5(f-fV 7 (x-x'). (3) 

In the context of the small-scale dynamo, this model was 
first proposed by Kazantsev (1967). While the Kazantsev- 
Kraichnan velocity field is, of course, highly artificial and does 
not approximate the real turbulent velocity field in any con- 
trolled sense, its performance in capturing the essential qual- 
itative and, in some cases, also quantitative, features of the pas- 
sive advection has been very impressive (on the passive scalar, 
see recent review by Falkovich, Gawedzki, & Vergassola 2001; 
on the kinematic dynamo, see Kinney et al. 2000, SCMM02, 
MCM02, Schekochihin et al. 2002d). This seems to suggest 
that the statistics of passive advection may be largely universal 
with respect to the structure of the ambient random flow. 

In the limit of large Pr m , the magnetic fluctuations are mostly 
excited deep in the subviscous range, in which the fluid mo- 
tions are strongly damped by viscous dissipation. The velocity 
field "seen" by the magnetic field is, therefore, regular and, in 
fact, effectively constitutes a single-scale flow. Most of the rel- 
evant statistical results turn out to depend just on the first few 
coefficients of the Taylor expansion of the velocity correlation 
tensor: 

K ij (y)=K Q 5'j- 1 - K2 (y 2 S i j-±y i y^ 

+\^y 2 [y 2 5 ii -\y i y^+--- (4) 

(some of the coefficients in the expansion are fixed by the in- 
compressibility constraint). Physically, k 2 ~ k v u is the stretch- 
ing rate, while ka sets the scale of the flow (m/^2 ~ k 2 ,) and is 
responsible for bending the field lines. 

2.2. Intermittency 

One of the simplest and most basic features of the kinematic 
diffusion-free regime is the lognormality of the distribution of 
the magnetic-field strength (see, e.g., Boldyrev & Schekochihin 
2000, SCMM02). On the most fundamental level, the lognor- 
mal character of the distribution of B follows already from the 
form of the induction equation (2) (without the diffusion term). 
Indeed, it is linear with respect to the magnetic field, which is 
multiplied by a random externally prescribed velocity-gradient 
matrix Vu. The Gaussianity of logB then follows by the Cen- 
tral Limit Theorem. 

Within the framework of the Kazantsev-Kraichnan model, 
this result can be derived exactly. The one-point probability- 
density function (PDF) of the magnetic-field strength in the 
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diffusion-free regime can easily be shown to satisfy the follow- 
ing Fokker-Planck equation: 



(5) 



Eq. (5) written in log variables is simply a ID diffusion equation 
with a drift and has the following Green's-function solution: 

-(1/2)« 2 I r°o dB , 



x exp 



[ln(B/B')-(l/4)K 2t y 

K2t 



(6) 



where Pq(B) is the initial distribution. The lognormality means 
that the PDF develops considerably spread-out tails at both 
large and small values of B, i.e., the fluctuating magnetic fields 
in the kinematic regime possess a high degree of intermittency. 
The moments (B n ) of the magnetic field have growth rates that 
increase quadratically with n: (B n ) oc exp[«(n + 3)K2f/4] . The 
exponential growth of the moments of B is the manifestation of 
the dynamo action of the turbulence. It measures the effect of 
stretching of the field lines by the ambient random flow. 

Note that these results apply to the diffusion-free situation 
(i.e., to an ideally conducting medium), and, therefore, hold 
during the time it takes the magnetic excitation to propagate 
from the velocity scales, where it is assumed to be initially con- 
centrated, to the resistive scales. This process occurs exponen- 
tially fast (Kulsrud & Anderson 1992; Schekochihin, Boldyrev, 
& Kulsrud 2002a), so the corresponding time can be estimated 
as t ~ Kj'logPr^ 2 , which amounts to several eddy-turnover 
times. The PDF of B with account taken of diffusion has not 
as yet been found, though Chertkov et al. (1999) did, under 
certain additional assumptions, derive the moments (B n ) in the 
diffusive regime. While the specific expressions for the growth 
rates of the moments are substantially modified, they still in- 
crease as n 2 , as would be the case for a lognormal distribution, 
so the intermittency is not diminished. 

2.3. The Folding Structure 

The geometry of the magnetic-field lines can be studied in 
terms of the statistics of their curvature K = b • Vb (here b = 
B/B). The evolution equation for the curvature can be derived 
from the induction equation: 



— K = K-(Vu) - (t-bb) 
at 



bKb: Vu-2Kbb: Vu 



+bb:(Wu)-(i-bb), (7) 

where we have dropped diffusion terms. For the Kazantsev- 
Kraichnan velocity, the PDF of curvature satisfies the following 
Fokker-Planck equation (SCMM02) 



, 7N d 1 20 

(l+K 2 ) P+ — P 

v ' dK K 7 



(8) 



where K has been rescaled by K* = (\2n\jl ni) x l 2 ~ k v . From 
this equation one immediately finds that the curvature attains a 
stationary distribution: 

(9) 



P(K)=l 



5 (1 + ^2)io/7- 

This PDF has a power tail ~ /f" 13 / 7 , which is in very good 
agreement with numerics (SCMM02). Naturally, if the resis- 
tive regularization is introduced, the power tail is cut off at 



K ~ k v ~ 'Px l J i 2 k v . The bulk of the curvature distribution is 
concentrated at values of curvature K <~ k v , which reflects the 
prevailing straightness of the field lines at subviscous scales. 
The field lines are significantly curved only in the bends of the 
folds. While these bends occupy a small fraction of the volume, 
there is a high degree of intermittency in the distrubution of the 
characteristic scales associated with them. This is reflected by 
the power tail of the curvature PDF. 

Furthermore, the field strength and the curvature are anti- 
correlated, i.e. the growing fields are mostly flat, while the 
curved fields in the bends of the folds remain relatively weak 
(cf. Drummond & Munch 1991; Brandenburg, Procaccia, & 
Segel 1995). It is not difficult to realize that the nature of 
this anticorrelation, which is derived as a statistical property 
in SCMM02, is, in fact, dynamical. Since this fact is impor- 
tant in our discussion of the nonlinear back reaction in § 3.2, 
we would like to explain it in more detail. 



Stretching of a fold 
by linear shear 




FIG. 2. — Stretching of a fold by linear shear. The bold arrows indicate 
the stretching direction of the shear. The dotted lines correspond to surfaces 
where B = 0. 

Let us examine what happens when a typical element of the 
folding structure (Fig. 2) is stretched by linear shear. Let Bf and 
Kf be the magnetic-field strength and curvature associated with 
the flat part of the fold, and B\, and K\, their counterparts in the 
curved part (the bend). In addition, let In be the length of the 
fold (parallel scale of the field), l± its thickness (perpendicu- 
lar scale of the field), and l\, the length of the bend. Fig. 2 can 
be thought of as a two-dimensional cross-section of a typical 
magnetic folding structure in three dimensions. It is assumed 
that the characteristic scale in the direction perpendicular to 
the page is approximately the same along the fold. Consider 
then the volume whose cross-section is indicated by the shaded 
area in Fig. 2: all the flux is through the surfaces whose cross- 
sections are represented by bold lines. Then, by flux conserva- 
tion, we must have 



Bfl± ~ ZVt>- 



(10) 



Suppose that the fold is stretched by a factor s in the parallel 
direction. Then lu —> sin, lb — > sib, an d, by volume and flux con- 
servation, Bf — ► sBf, so the field is amplified in the flat part of 
the fold. We now observe that Kf~l/ Z|| and Kb ~ whence 
B f K f ~ B f //|| and, using Eq. (10), B b K b ~ B h /l± ~ B f /l b . Both 
of these ratios remain unchanged during the stretching event. 
We conclude that linear shear transformation preserves 



KB ~ const 



(11) 
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FIG. 3. — Instantaneous cross-sections of the field strength B (left) and curvature K (right) during the nonlinear stage of a simulation with v = 5 ■ 10 2 , r) = 2 ■ 10 
(Pr m = 250). Darker areas correspond to larger values. Detailed anticorrelation between B and K is manifest. 



throughout the fold: hence, the anticorrelation between B 
and K. Note that Eq. (11) can be regarded as equivalent to the 
conjecture by Brandenburg, Procaccia, & Segel (1995) that flux 
tubes straighten at the same rate as the corresponding magnetic 
fields grow. 

The perceptive reader might wonder about the consistency 
of Eq. (11) with our previously stated results that, while the 
statistics of the field strength are lognormal, the curvature dis- 
tribution has a power tail. The explanation is as follows. The 
argument that led to Eq. (11) was based on approximating the 
ambient flow by a linear shear and thus described the effect of 
the first derivatives of the velocity on the magnetic field. We 
now observe that the curvature evolution equation (7) also con- 
tains the second derivatives of velocity, which enter as a source 
term. Their role is to bend the field lines at the scale of the flow 
(see Fig. 1): once the field becomes too flat (i.e., flatter than 
the flow), it is bent within one eddy-turnover time and K <~ k v 
is restored. Therefore, the curvature cannot remain below val- 
ues of order k v for a long time (or in many places). In formal 
terms, the second derivatives of velocity break the scale invari- 
ance of Eq. (7) and, consequently, of the associated Fokker- 
Planck equation (8) (this is reflected in the presence of unity in 
the first term on the right-hand side). A stationary distribution 
with a power tail is thereby made possible. In fact, this power 
tail is a limiting envelope for a lognormal tail produced by the 
action of the linear shear (see SCMM02). 

The most important consequence is that the field settles into 
a reduced-tension state: the tension force can be estimated by 
B VB - k\\B 2 ~ KB 1 ~ k v B 2 on the average. A simple reductio 
ad absurdum argument can be envisioned to further support this 
statement. Let us write the evolution equation for F = B • VB: 



indeed, 



dt 



F = F Vu + BB : VVu, 



(12) 



where the diffusion terms are again dropped. Suppose for 
a moment that the field is chaotically tangled, i.e. F ~ k\\B 2 
with fey ~ti> k v . Then the term in Eq. (12) that involves 
the second derivatives of the velocity field can be neglected and 
Eq. (12) becomes formally identical to the evolution equation 
for B. The moments of F must, therefore, grow at the same rates 
as the moments of B, and we estimate (F 2 ) /(B 4 ) oc (B 2 ) /(B 4 ), 
which decays exponentially fast in time. An exact statistical 
calculation for the Kazantsev-Kraichnan velocity, shows that, 



(|B-VB| 2 ) 

(B 4 ) 



28 K4 

~9^2 



k 2 



(13) 



asymptotically with time, starting from any initial conditions 
(SCMM02). The convergence is exponentially fast at the 
stretching (eddy-turnover) rate. We conclude that even an ini- 
tially chaotically tangled magnetic field will quickly develop 
the folding structure. 

Thus, the nonlinear saturation, which is due to the Lorentz 
tension balancing the stretching action of the flow, occurs when 
the energy of the field becomes comparable to the energy of 
the turbulent eddies. Note that in a hypothetical chaotically 
tangled field with h\ <~ k±, the tension would be much larger: 
B • VB ~ k ri B 2 , so saturation would be possible already at very 
low magnetic energies. 

3. THE NONLINEAR REGIME 

3.1. Numerical Results 

No satisfactory analytical description of the nonlinear state 
is as yet available, so one must be guided by results of numer- 
ical experiments. The main obstacle in the way of a defini- 
tive numerical study is the tremendously wide range of scales 
that must be resolved in order to adequately simulate the large- 
Pr MHD: indeed, one must resolve two scaling intervals, the 
hydrodynamic inertial range and the sub viscous magnetic one. 
Since this is not feasible, we propose to simulate the initial 
stage of the nonlinear evolution up to the point when the to- 
tal energy of the magnetic field equalizes with the energy of 
the viscous-scale turbulent eddies. This stage can be studied 
in the viscosity-dominated MHD regime where the hydrody- 
namic Reynolds number Re is order one and the external forc- 
ing models the energy supply from the larger eddies (cf. Cat- 
taneo, Hughes, & Kim 1996; Kinney et al. 2000). Moreover, 
one can argue (MCM02, Schekochihin et al. 2002b) that once 
the magnetic energy does equalize with that of the smallest ed- 
dies, the following scenario takes place. The magnetic back 
reaction leads to suppression of the shearing motions associ- 
ated with the viscous-scale eddies. Larger-scale eddies, which 
are more energetic (but have slower turnover rates), continue to 
drive the small-scale magnetic fluctuations by the same stretch- 
ing mechanism that the viscous-scale ones did, until the mag- 
netic field becomes strong enough to suppress these eddies as 
well. This process continues until all field-stretching motions 
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throughout the inertial range are suppressed (see Schekochihin 
et al. 2002b, for further discussion). Speculatively, this could be 
thought of as some effective renormalization of Re, so that the 
final statistics of the magnetic field would again be described by 
the low-Re MHD model. Some numerical evidence supporting 
this picture is given by MCM02. 

In our simulations, we choose the parameters so that the forc- 
ing and the viscous scales are comparable and the statistics of 
the subviscous-range magnetic fluctuations can be studied al- 
ready at 128 3 resolution. The numerical set-up and the spectral 
MHD code used are exhaustively described in MCM02. The 
forcing is large scale, nonhelical, and white in time. The units 
are based on the box size 1 and the forcing power 1 . In these 
units, setting v = 5 • 10~ 2 effectively leads to Re ~ 1. 

After the initial kinematic growth stage, saturation is 
achieved where the total magnetic and hydrodynamic energies 
are equal. 1 We then measure the PDFs of the magnetic-field 
strength and curvature in the saturated state. Our findings are 
as follows. 



-_ i | i | i | i | i | i | i | i | i r 
: PDFs of field strength normalized to B 




FIG . 4. — PDFs of the field strength: onset of nonlinearity and saturation. The 
field strength is normalized to its rms value for each PDF. The intermittency at 
earlier times is thus clearly seen to be stronger than during the saturated stage. 



1 . The level of intermittency is lower than in the kinematic 
case, the field-strength distribution developing an exponential 

tail. The kurtosis of the field decreases to (B 4 )/(B 2 ) ~ 4 after 
growing exponentially during the kinematic stage (for a Gaus- 
sian three-dimensional random field, the kurtosis would have 
been 5/3). 

The partial suppression of intermittency indicates a plausible 
scenario for the onset of nonlinearity (see Fig. 4). As we saw 
in § 2.2, the kinematic dynamo gave rise to a highly intermit- 
tent lognormal spatial distribution of the field strength. Already 
early on, when the total magnetic energy is only a fraction of 
its saturated value, there are tiny regions where the magnetic- 
energy density locally approaches values comparable to, and 



greater than, the energy density of the fluid motions. This acti- 
vates the nonlinear back reaction in these regions, so as to sup- 
press further growth of the field there. As the overall amount 
of the magnetic energy grows, the fraction of the volume where 
the nonlinearity is at work increases until the nonlinear state 
is established globally. In the process, the magnetic fields be- 
come more volume-filling and, consequently, less intermittent 
(cf. Brandenburg et al. 1996; Cattaneo 1999a,b). The final PDF 
in the fully nonlinear regime has an exponential tail (Fig. 4): 
still considerably intermittent, but less so than the lognormal 
distribution. 
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FIG. 5. — PDFs of the field-line curvature in the nonlinear regime. The solid 
curve is for a simulation that used 8th-order hyperdiffusion. At curvatures ap- 
proaching resistive cutoff, K ~ k v , the scaling is destroyed. 

2. The folding structure of the field is unchanged from the 
kinematic case. The anticorrelation between the field strength 
and the curvature persists (see Fig. 3, which illustrates the dy- 
namical nature of this anticorrelation discussed in § 2.3). Their 
correlation coefficient r Kfi = (K 2 B 2 ) / (K 2 ) (B 2 ) - 1 ~ -0.7 in our 
simulations, which is quite close to its minimum allowed value 
of — 1. Even more remarkably, the distribution of the curvature 
retains the same power tail ~ K~ 13 / 7 (Fig. 5). We find this last 
feature particularly striking: indeed, not only does a kinematic 
theory based on the synthetic Kraichnan velocity field predict 
a quantitatively correct nontrivial scaling for the tail of the cur- 
vature PDF, but this scaling also survives in the fully nonlinear 
case! 

Finally, we would like to offer the following caveat with re- 
gard to the numerical study of the small-scale magnetic turbu- 
lence by means of spectral methods. We have found that if the 
amount of dissipation applied to the magnetic field is not suf- 
ficiently large to completely override the effect of dealiasing 
(see Canuto et al. 1988), the folding structure is destroyed and 
incorrect curvature statistics are obtained (cf. MCM02). The 
resulting resolution constraints tend to be quite stringent: it is 



See the argument at the end of § 2.3. The equipartition only holds asymptotically with Pr m . When the magnetic diffusivity is too large, saturation occurs at smaller 
magnetic energies (for more details, see Schekochihin et al. 2()02b,d). This probably accounts for subequipartition saturated states seen in some of the simulations of 
small-scale dynamo (e.g., Brummell, Cattaneo, & Tobias 2001). We note, however, that the lower magnetic-energy levels in the saturated state do not alter either the 
intermittency- or the structure-related properties of the small-scale fields. 
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not enough to simply make sure that the UV cutoff of magnetic- 
energy spectrum is resolved. We believe that the destruction of 
the structural properties of the field by dealiasing is due to non- 
local nature of the corresponding operator. Consistently with 
this view, hyperdiffusion does not exhibit the same adverse ef- 
fect (see Fig. 5). 

3.2. A Model of Back Reaction 

The physical reason for the preservation of the folding struc- 
ture is, of course, that the fluid motions at subviscous scales are 
strongly damped and cannot "unwrap" the folds. Here we use 
the folding picture to devise a simple ad hoc physical model that 
explains the two key quantitative results that have emerged from 
our numerics: the exponential tail of the PDF of the magnetic- 
field strength and the unchanged kinematic power tail of the 
curvature PDF. 

It is clear that, once the magnetic fields grow strong enough, 
they will resist further stretching by the ambient velocity shear. 
This increased rigidity of the field can be modelled by adding a 
nonlinear relaxation term to the induction equation: 

^B = B-Vu ext -r; 1 (B)B, (14) 
dt 

where u ex t is the part of the velocity field due to the external 
forcing and 7>(B) is some effective nonlinear relaxation time. 2 
In view of the heuristic way in which the nonlinearity has 
been introduced into the induction equation, Eq. (14) cannot be 
expected to, and, indeed, does not, correctly model the spatial 
properties of the magnetic field. What it effectively amounts to 
is a simple local one-point closure. If the external velocity field 
is modelled by the Kazantsev-Kraichnan velocity, the Fokker- 
Planck equation (5) derived in the kinematic case can easily be 
amended to include the nonlinearity: 



1 



■ K2-S 



d 1 



p+t;\b)p 



4 ' "" dB B 2 ' 
This equation has one normalizable stationary solution 



P(B) = const B 2 exp 
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dB 1 



(15) 



(16) 



where B t y P icai ~ B ms is the typical value of the back-reacting 
magnetic field in the system. This regime corresponds to what 
is sometimes referred to as the the viscous relaxation effect 
(Chandran 1998; Kinney et al. 2000). Substituting the expres- 
sion (19) for t~ 1 (B) into Eq. (16), we find a Gaussian PDF (cf. 
Boldyrev 2001) with B rms ~ {k,2v) 1 I 2 , which approximately cor- 
responds to magnetic energy equalizing with the energy of the 
viscous-scale eddies and accounts for the energy equipartition 
we have seen in our simulations. 

However, these arguments are only appropriate in applica- 
tion to the statistics of the typical values of the magnetic- 
field strength, i.e., to the core of the distribution function (cf. 
Boldyrev 2001). In order to understand the tail of the PDF, one 
must consider rare events, specifically, the instances (or places) 
where B is atypically large. In the framework of the "stretch- 
and-fold" small-scale dynamo mechanism, such events can oc- 
cur if magnetic field gets "overstretched" (for instance, because 
the shear has acted coherently in one direction for an atypically 
long time). The local relation between the curvature and the 
field strength associated with it is given by Eq. (11), whence 
KB ~ fc„fitypicai- Inserting this into Eq. (18), we find 

t _1 {B) _ ^typical ^ B ^ B ^ ^ (2Q) 



K2Jo B' 

In order to make an intelligent guess about the B-dependence 
of t~ 1 (B), we have to invoke our understanding of the field 
structure at subviscous scales. The relaxation term in Eq. ( 14) is 
due to the suppression of the velocity shear by the Lorentz back 
reaction. The latter can be thought of as inducing a certain mag- 
netic component u m of the velocity field that counterbalances 
the external shear, so t~ 1 {B) <~ Vu m . Since we are consider- 
ing a viscosity-dominated regime, Vu m should be estimated by 
balancing the viscous damping and the magnetic tension: 

-j/Au m -B-VB. (17) 
The characteristic scale of u m cannot be smaller than the vis- 
cous scale: otherwise, u m would be viscously damped at a 
shorter time scale than that of the external shear and could not 
provide an effective counteraction to it. We have, therefore, 

, B VB KB 1 

T;\B)~k v ii m — . (18) 

vk v vk v 

We emphasize that the above relation represents a local force 
balance in a particular fold. Now, in a typical case, the folding 
structure would imply K ~ /Cjy, SO 

B 2 

T;\B)~—, B ~ ^typical, (19) 

2 A model of nonlinear feedback mathematically very similar to this one was studied 
field in a random medium. 



which, upon substitution into Eq. (16), gives the exponential 
tail evidenced by the numerics. 

The persistence of the kinematic curvature statistics in the 
face of nonlinear effects can also be understood in these terms. 
Essentially, the tail of the curvature PDF remains unaffected 
by the back reaction because it describes areas of anomalously 
large curvature (K^k^) where, due to the anticorrelation prop- 
erty, the field is weak. On a slightly more quantitative level, 
we argue that the effect of the back reaction on the curvature 
can also be modelled by a simple nonlinear relaxation term, as 
inEq. (14): 

— K= [rhs of Eq. (7) with Vu e xtl -t;\B)K. (21) 
dt 

Here t~ 1 {B) is again estimated via Eqs. (18) and (1 1). The non- 
linear relaxation term in Eq. (21) is then 



-l, mr {KB) 2 „ ^^typical ~ 

-t/(B)K — n ^^n, 

vk v v 



(22) 



where n = K/K. With this correction, the Fokker-Planck equa- 
tion (8) for the PDF of curvature becomes 



7 d 
d t p=- Kl —K 
4 dK 



d 1 



20 



l+K z ) P+—P--P 



dK K 



K 



, (23) 



where a is a numerical constant of order unity and, as before, 
K is rescaled by K* <~ k v . It is a straightforward exercise to 
show that the stationary PDF is now given by 



P(K) = const 



]£ e a arctan(/Q 



(l + K 2 ) 



10/7 : 



(24) 



which has the same power tail <~ K~ n l 1 as its kinematic coun- 
terpart (9). 

Finally, we would like to stress the qualitative character of 
the ideas and results put forward in this section. A more quan- 
titative nonlinear theory based on these ideas may be feasible, 

by Zeldovich et al. (1987) in the problem of self-excitation of a nonlinear scalar 
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but is left for future work. Another important issue that requires 
careful quantitative treatment is the role of Ohmic diffusion. 
Based on the model of back reaction proposed here, we seem 
to be able to understand the nonlinear regime without including 
the resistive terms. These terms are hard to treat analytically 
due to the usual closure problem associated with the diffusion 
operator. Numerically, we have confirmed that the field struc- 
ture is unaffected by a switch to hyperdiffusion (see Fig. 5), 
but the locality of the dissipation operator probably does matter 
(see remarks at the end of § 3.1). Further investigation of the 
universality of the statistics of the small-scale magnetic turbu- 
lence with respect to the form of the UV regularization is also 
left for the future. 

4. CONCLUSIONS 

We have found that the folding structure of subviscous-scale 
magnetic fluctuations that is formed via the kinematic "stretch- 
and-fold" small-scale-dynamo mechanism remains the essen- 
tial feature of the nonlinear regime. The scale separation is of 
crucial importance here: while small-scale structure can be gen- 
erated and maintained in large-scale random flows, these flows 
lack the ability to coherently undo the structure even when act- 
ing in concert with the magnetic back reaction. 

Both our theoretical arguments and our numerical experi- 
ments were based on viewing the turbulent velocity field as a 
single-scale random flow. In real turbulence, Re is, of course, 
fairly large (10 4 in the ISM), so many hydrodynamic scales 
come into play. Unfortunately, the resulting scale ranges are 
too broad to be adequately simulated. It is interesting, however, 
that the results presented above appear to hold in simulations 
with more realistic Re (up to 10 3 ) but relatively small Pr (down 
to values of order one), in accordance with the arguments pre- 
sented at the beginning of Sec. 3. We believe, therefore, that 
the physics and the numerics we have laid out provide at least 
a semiquantitative description of the statistical properties of the 
small-scale MHD turbulence in high-Pr m plasmas. 

The implications for astrophysical objects can be significant. 
For the large-scale galactic dynamo, small-scale fields must be 



taken into consideration if a nonlinear aft theory is to be con- 
structed. If the net effect of the accumulated small-scale mag- 
netic energy is to suppress a, an alternative theory will have 
to be sought. The nonlinear evolution of the small-scale mag- 
netic turbulence in protogalaxies (Kulsrud et al. 1997) and in 
the early Universe (see, e.g., Son 1999; Christensson, Hind- 
marsh, & Brandenburg 2001) determines the energy and the co- 
herence scale of the seed field inherited by newly formed galax- 
ies. The structure of tangled magnetic fields in the intracluster 
gas crucially affects thermal conduction in the galaxy clusters 
(see, e.g., Chandran & Cowley 1998; Malyshkin 2001; Narayan 
& Medvedev 2001). In the accretion-disk physics, presence of 
large amounts of small-scale magnetic energy could lead to new 
models for the angular-momentum transport (Balbus & Haw- 
ley 1998) and for the acceleration of jets (see, e.g., Heinz & 
Begelman 2000). In the solar astrophysics, the possibility was 
recently raised that a substantial part of the magnetic energy in 
the quiet photosphere of the Sun resides in small-scale mag- 
netic fluctuations (Cattaneo 1999b). Indeed, it is natural to ex- 
pect that, just as turbulence itself, small-scale random magnetic 
fields are ubiquitous in the Universe. 

Finally, constructing a self-consistent physical theory of the 
small-scale magnetic turbulence constitutes a fascinating task 
in its own right. Though fifty years have passed since Batch- 
elor (1950) took the first steps down this road, an inquisitive 
researcher will still find surprises at every turn, and it might 
well be short-sighted to claim that we are able to discern the 
contours of the final destination. 
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Grants No. AST 97-13241 and AST 00-98670, and the USDOE 
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